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Abstract 

Biharmonic curves are a generalization of geodesies, with apphcations 

in elasticity theory and various branches of computer science. The pa- 

^J I per proposes a first study of biharmonic curves in spaces with Finslerian 

geometry, covering the following topics: a deduction of their equations, 
existence of non-geodesic biharmonic curves for some classes of Finsler 
2 . spaces and specific properties. We prove, for instance, that geodesic cur- 

vature is constant along Finslerian biharmonic curves. Integration of the 
biharmonic equation is presented for two concrete Finsler metrics. 
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■^ . 1 Introduction 

■"sj" , In elasticity theory, as well as in other domains, such as computer graphics and 

^^ ' image processing, a classical problem is that of finding the extremal curves of 

the elastic energy (ot bending energy) functional, defined as the integral of the 
squared curvature k^ of a curve with respect to its arc length s. Its extremal 
curves (called least- energy or elastic curves) are interpreted as the "smoothest" 

k> ' curves joining some given points or with given boundary conditions, [1], [5], [8]. 

}_( , The Riemannian generalization of the elastic energy, called the bienergy, is 

.^.' defined as: 

E2ic) = J \\T{c)f ds, (1) 

c 

where t{c) is the tension of the curve, [H], [2]. In the natural parametrization, 
we have t{c) = VcC - that is, the integrand is actually the squared geodesic 
curvature k^ = l|Vcc|| - and extremal curves of the bienergy, i.e., biharmonic 
curves, are described by the equation: 



V^c = R{c, ^cc)c. (2) 

Riemannian spaces, bih 
maps 



In Riemannian spaces, biharmonic curves - and, more generally, biharmonic 
ips - were examined quite in detail, |3, [6], [7], [10], [12], [13], [M]. We will 



only mention here several properties. First of all, geodesies of V are always bi- 
harmonic (actually, they are minimum points for the bienergy) , but the converse 
is generally not true. 

A non-harmonic biharmonic curve is called proper biharmonic. The problem 
of the existence of proper biharmonic curves is one of the most discussed in the 
literature. Thus, f7\. [T^: 

- In the Euclidean space R" , every biharmonic curve is a segment of a straight 
line. 

- On surfaces with nonpositive Gauss curvature, any biharmonic curve is a 
geodesic. 

- In Riemannian spaces with nonpositive sectional curvature, every closed 
biharmonic curve is a geodesic. 

Also, along any (unit speed) biharmonic curve, the geodesic curvature k is 
a constant. 

Passing to Finsler spaces, the interest in biharmonicity is motivated, beyond 
the purely mathematical reasons (a list of which can be found, for instance, in 
[H]), by the potential of Finsler geometry for models in elasticity theory, [2], 
[15] . In the present paper, we find the equations of biharmonic curves in Finsler 
spaces and study the Finslerian counterparts of the above mentioned results. As 
a somehow surprising fact, geodesic curvature remains constant along unit speed 
biharmonic curves. Still, to a difference from the Riemannian case, in dimension 
greater or equal to three, there exist flat (locally Minkowskian) Finsler spaces 
which admit proper unit speed biharmonic curves. 

In the last section, we determine proper biharmonic curves for two concrete 
Finsler metrics. 

2 Finsler structures 

In this section, we present known notions and results in Finsler geometry to 
be used in the following. Consider a C°°-smooth, n-dimensional differentiable 
manifold M and denote by {x,y) :— {x^,y^)i^Yn^ ^^^ local coordinates on its 
tangent bundle {TM,Tr,M). By ^i and .i we will mean partial differentiation 
with respect a:* and y' respectively; by r{E) we will mean the module of the 
sections of some bundle E over M. 

A Finsler structure, W, f^, on the manifold M is a function F : TM — s> R 
with the properties: 

1) F{x,y) is C°°-smooth for y y^ and continuous at y = 0. 

2)Fix,Xy)^XF{x,y),yX>0. 

3) The Finslerian metric tensor: 

g^j{x,y):^ -{F'^{x,y)).ij (3) 

is positive definite. 

By its dependence on the fiber coordinate j/, the Finslerian metric tensor gij 
"lives" on the tangent bundle TAf, hence it is natural to perform the study on 



this bundlcEl, [2], El- 

The arc length of a curve c on the Fmsler space {M,g) is given, [3], by: 

l{c)= f F{x,dx). (4) 



A key ingredient in Finsler geometry, which considerably simplifies cal- 
culations, are adapted frames given by Ehresmann (nonlinear, [2], [9], |4]) 
connections on TM. The typical choice is the Cartan nonlinear connection 
TTM = HTM ® VTM, having the property that complete lifts to TM of 
tangent vectors to unit speed geodesies of M are always horizontal. The Cartan 
nonlinear connection is locally defined by the coefficients 

G j '■— G .J, (5) 

where G* — G'(x, y) are obtained from the equations of unit speed geodesies of 

{M,F): 

^ + 2G^(a;,y)=0, y ^ x. (6) 

The (generally, nonholonomic) adapted basis on T{TTM) is then given by: 

(,.._!, -GVx,.)|j. ft^A) (T) 

and its dual is: (da;*, 5y^ — dy^ + G^^dx^). Any vector field X on TM can 

be decomposed as: X = hX + vX, hX — X^Si, vX = X^di\ its horizontal 
component hX and its vertical component vX are vector fields on TM . 

In the presence of a nonlinear connection, vector fields on M can be identified 
with horizontal vector fields on TM via the horizontal lift 

Z'' : r(rAf ) -^ T{HTM), X = X.'d, ^ l''{X.) = (X' o 7r)5,. 

The Finslerian metric tensor g gives rise to a scalar product 

{X, Y) = g,,X'Y^ 

of horizontal vector fields X = X'd^, Y = Y^Sj e T{HTM). 
The following (1,1) -type tensor fields on TM: 

J = di® dx\ e = 5i® 5y\ (8) 

have the role of mapping horizontal vector fields into vertical ones and vice- versa 
{J{5,) = 5„ J{di) = 0; e{d,) = 5„ e{5^) = 0). 



^Another natural variant is the pullback bundle n*TM, used in |11) . [2]. 



As covariant differentiation rule for tensors on TM^ we will use the one 
given by the Cartan affine connection D on TM, [2], [9], locally described by 
the coefficients: 

Dsj, = r,J^, DsA - r'^.,a„ d^^6, = c^.-^^ ^4 4- = C^3k^^, (9) 

where T*^^ = V ^f.{x,y), C*^^ = C*j^(a;, y) are as follows: 

The Cartan connection is metrical and preserves by parallelism the horizontal 
and vertical distributions on TTM , i.e., 

DxihY)^hDxY, DxivY) = vDxY, yX,Y er{TTM). (11) 

Its torsion T has, in the adapted frame, the following nonzero local compo- 
nents: 

r(X, Y) = W^^d^ ® dx^ ® dx' + C'jk^, » Sy'' (g) dx^ + P^^d^ (g) Sy'' (g dx^ , (12) 

with C^jj, as in (fTO)) and: 

R^^,^6,G^^-S,G\, P),=y'C^^,^^ (13) 

(where |; := -D^J. The first of the three terms in p^ is the curvature of the 
nonlinear connection N and the other two define two specific geometric objects 
to Finsler geometry: 

2) the CaHan tensoB- C{X,Y) := hT{JX,Y), yX,Y e T{HTM), in local 
writing: C = d^jk^i ® dx'' (g dx^ ; 

3) the Landsherg tensor P{X,Y) := er{JX,Y), yX,Y e T{HTM), in 
local coordinates: P := P^jk^i ^ dx'' ® dxK 

The Cartan and the Landsberg tensors obey similar properties: 
i) the scalar products (C(X, F), W) , {P{X, Y), W) are totally symmetric in 
X,Y,W eT{HTM); 

ii) contraction of C and P with y provides zero, i.e.: 

C(/"(y), 0=0, P(/"(2/), 0=0. (14) 



The curvature TZ of the Cartan connection D also splits in the adapted 
basis into horizontal and vertical components. The local expressions of these 
components can be found, for instance, in [5], [5]; here, we will only mention 
some of their properties: 

a) VC/, V G T{TM) : Tl{vU, hV){l'^y) = P{eU, hV); TZ{vU, vV){l'^y) = 0; 



■^The Cartan tensor is a measure of the non-Riemannian character of the Finsler metric; 
the equality C = is equivalent to the fact that the metric is Riemannian. 



b) the h{hh) component of TZ gives rise to the Jacobi endomorphism of 
{M,F), [2], understood here as a (l,l)-type horizontal tensor: 

R{X) = n{X,l''y)l''y, VX £ T{HTM); (15) 

in local writing, R = R^ j5i ® dx^ , i?' = i?'^ j/'^. The components Rik are 
symmetric, therefore, 

{R{X),Y) = {R{Y), X) , VX, y e T{HTM). (16) 

If the system {l''{y),X} is orthonormal, then K{X) := —{R{X),X) is 

equal to the flag curvature of (Af, i^), [4], which represents an analogue of sec- 
tional curvature from Riemannian spaces. 
Particular cases: 

1) In Riemannian spaces, we have gij — gij{x), C — 0, P — 0, while the 
h{hh) component n{hU, hV)hW {U, V,W e r{TM)) is, up to a horizontal lift, 
the Riemannian curvature of {M,g). 

2) Locally Minkowski spaces ("flat" Finsler spaces) are spaces in which there 
exists around any point a local chart in which gij = gij(y). In this case, F* j, = 0, 
R = 0, P = 0, but, generally, C 7^ 0. Geodesies of a locally Minkowski space 
are straight lines. 

3) Landsberg spaces are defined by the condition P — 0. 



3 Bienergy of a curve and its first variation 

Consider a curve c : / ^' M, s 1— >■ (a;*(s)) (where / C M is an interval), 
parametrized by its arc length and its lift c' : I —^ TM, s n- {x^ (s) , x^ (s)) 
to TM; that is, c has unit Finslerian speed: 

F{x,x)^l. (17) 

We denote simply by D :— D'^ the puUback of the Cartan afBne connection D 
bye'. 

Tangent vectors of c' are expressed in the adapted basis as: 

T':=c'M) = y'5, + "-j-d,, y' ^ x\ 
as 

Taking into account that F* j.j/'' = G" ([2], [3]), the components of T' in 
the adapted basis can be re-expressed as: 

T := hT' = l''{y), vT' = JiDgT) (18) 

(where covariant differentiation here and in the following is made, [4], with 
"reference vector" T, i.e., in the expressions of the connection coefficients F = 

r{x, y) and C — C{x, y), we set y' = i*). 



The energy of the curve c is given, |llj , by 

Ei{c) = y {T,T)ds. (19) 

c 

Extremals of the energy, i.e., harmonic curves, |11) . are constant speed geodesies 
of (M, F), described by: 

Do.T^O. (20) 

The horizontal vector field t{c) :— DqT, defined along the curve c, is called the 
tension of c. Thus, it appears as natural: 

Definition 1 1) The bienergy of a curve c on a Finsler space (M, F) is: 

E2{c) = iy" {DqT, DgT) ds. (21) 

c 

2) Critical points of the bienergy are called biharmonic curves of{M,F). 

There holds the following result: 

Theorem 2 A curve c on the Finsler space {M, F) is biharmonic if and only 
if: 

T2{c) := DlT + RiDeT) - P{DoT, DgT) + C{DgT, DIJ) = 0. (22) 

Proof. Let us consider smooth variations Ce{s) = c[e,s), cq{s) = c{s) with 
fixed endpoints of c and their lifts c^ = (c^, — ^) to TM. With the notations: 

T' := 5t(a,), V := c',{d,)- T := hT , V := W, (23) 

we have vT' = J{DaT) and, at e = 0: 

T = Tje=o, V^:=VJe=o. 
Taking into account that the Cartan connection D is metrical, we can write: 

^{c,) = ^J-^{D9T,DgT)ds^ J {Do^DaT,DoT)ds. (24) 

c c 

Further, since D preserves the distributions generated by A^, (1^51) . 

Dg^DoT = n{V',T')T+ Da^Dg^ihT) =n{Y',T')T+ 
+ DoShDoX) -7^(V',T')T+ DQ^hTiV ,T)) + Dl^{h\')- 

with V = KV' and substituting into (|24)) . we get: 

^ = / (7^(V^ T')T + i^a, (/ir(V', T')) + Dly, DgT) ds. (25) 



Let us evaluate each term in the above integraL First of ah, the curvature term 
is: 

/i := /'(7^(V',T')T, DoT)ds= f {n{V,T)T + n{vV' ,T)T + 

c c 

+ n{V,vT')T + n{vV',vT')T, D9T)ds. 

Along each of the curves c^, we have T = l^{y)\ taking into account the 
properties of the curvature components, the terms (7?,(uV', T)T, Dg^T) and 
((R,{vV\ uT')T, Dg^T) vanish, while the two remaining ones can be written in 
a simpler form; 

{TZ{V,T)T,D9T) - {R{V),D9T)^ {R{DsT),V); 
{TZ{hV',vT')T,D9T) = -{P{D9T,V),D9T))^~{P{D9T,D9T),V)) 

therefore, we can write: 

h= f (RiDaT) - P{D9T,D9 T),Y) ds. (26) 

C 

The second (torsion) term in (PSJ) . i.e.: 

h- f {D9Shr{V',T')), D9T)ds 



will be integrated by parts: 

h - J{ds {hr(V',T'),D9T) - ((/.r(V',T')), DlT)}ds. 

C 

Since the variation has fixed endpoints, we have: Jds {hT(V\ T'), Dq/T) ds - 

c 

0; using (fT^ and (fH)) in the evaluation of the remaining term, we are led to: 

h= j {C{D9T,Y), DlT) ds = J {C{D9T,DlT), V) ds. 



Finally, integrating twice by parts the third term in (|25p . i.e., 

h:= J {DlY, D9T)ds, 

c 

we are led to: 



h:^ J{DlT,V)ds. 



Summing up and evaluating at e = 0, the first variation of the bienergy is: 

de 
= / {DlT + R{DoT) - P{DaT, DqT) + C{DsT, DIT),V) ds, 

c 

which proves the statement. ■ 

The vector field T2(c) in ((22)) is called the bitension of c. In the following, 
we will also write the biharmonic equation (|22p as: 

D^T^A, (27) 

where 

A := -R{DqT) + P{DqT,DqT) - C{DqT,DIT). (28) 



Remark. Any unit speed geodesic on M is biharmonic; a non-geodesic 
biharmonic curve will be called proper biharmonic. 

Particular cases: 

1) If M is a Riemannian space, then the equation of biharmonic curves 
reduces to: 

Dl^ T + R{Do, T) = ^ Dl^ T-n{T,DaT)T = 0. 

Up to a projection onto the base manifold, this is the known equation ^. 

2) In locally Minkowski spaces, equation (|22|) becomes: 



DlT + C'iDgT, DlT) = 0. 
3) In Landsberg spaces, T2(c) = D^T + RiDgT) + C{DqT, DIT). 



4 Frenet frame and geodesic curvature 

We will suppose that, along c, the vectors {T, DqT, Dg T, ..., Dg T} are linearly 
independent; in particular, this implies that Dg^T ^ 0, i.e., c is non-geodesic. 
Similarly to Riemannian geometry, [3] , |10) , we get the Frenet frame along c as 
the orthonormalization of these vectors, obeying the relations: 

f ei = T = l'\y) 
Dd.ei = Kie2 

^9,62 = -Kiei + K2e3 , . 

Da.Si^ '~Ki^iei^i + KiCi+i, i = 2,..,n-l 



here, ei,...,e„ are regarded as sections of the bundle (c') ^{HTM)); the first 
curvaturqfl ki, i.e., 



Ki = WDqTW := ^{DaT,DgT) (30) 

is called the geodesic curvature of c. 

Using (P5|) . the expressions of the covariant derivatives Dg T, i = 1, 3 in the 
Frenet frame are: 

Dd^T — Kie2 

Dg^T — — KiCl + K[e2 + KlK2e3 

Z?| r = —3KlK[ei + (k" — K^ — KlK\)e2 + (2k'j^K2 + KlKjjCS + KlK2'«3e4. 

(31) 
Substituting (pij) into the biharmonic equation (P7| and using the fact that 
the Cartan tensor C vanishes on ei, we have: 

A = -^1^(62) + KiP(e2, 62) - KlK'iC(e2, 62) - K^K2C(e2, 63). (32) 



Taking the scalar product of the third equality pip with ei,i — 1, n and noticing 
that {A, ei) =0, we find: 

Proposition 3 A unit speed curve c : I ^ M is biharmonic if and only if: 

kik[ = 0; 

k'( — k\ — kik\ = {A, 62) 

2k'iK2 + ki4 = (^1,63) (33) 

K1K2K3 = (^,64) ; 
^ 0= (^e^), i> 5. 



The assumption that c is non-geodesic means that ki 7^ 0, which, by (|33]), 
gives us k'i = 0. That is, if c is a proper biharmonic curve parametrized by arc 
length then: 

Ki — const. 7^ 0. 

Since in the case when ki = (i.e., c is a geodesic), the curve is trivially 
biharmonic, we get a similar result to the one in Riemannian geometry: 

Theorem 4 Along any biharmonic curve of a Finsler space, the geodesic cur- 
vature Ki — \\Dg^T\\ is constant. 

5 Existence of proper biharmonic curves 

5.1 Closed biharmonic curves 

We will assume in the following that c is a smooth closed biharmonic curve, i.e., 
c and its derivatives coincide at its endpoints. We can therefore write: 

= J^^ {D9T,DaT)ds = j {{D%T,DqT) + {DlT,DlT)}ds. 



'in dimension 2, the irst curvature is replaced by the signed curvature, [12) 



Using the biharmonic equation, this is: 

0= J{{A,DsT) + {DlT,DlT)}ds. (34) 

C 

Let us evaluate J {{A, Dq^T) ds. By virtue of the total symmetry of C, the 

c 

term involving it in the mentioned integral is: 

j{-C{DaT,DlT),DeT)ds = -j{C{DeT,DoT),DlT)ds. (35) 

c c 

Since the curve is closed, we have jds {{C{Dq^T,DqT)),DqT) — 0, hence, 

c 

integration by parts in ([55)) leads to: 

/ {-C{D9T,DgT),DlT) ds = J {DoAC{D9T,DgT)),D9T)ds = 
' = / {{Da^QiDgT, DqT), DqT) + 2 {C{DaT, Dlj), D9T)}ds 

c 

Consider the following tensor acting on horizontal vector fields: 

C{X,Y,Z):={DjxC){Y,Z), yX,Y, Z eT{HTM); (36) 

then, using the total symmetry of C and (J13p . we can write: 

/ {-CiDgT, Dlj), DqT) ds = \j {{Do^OiDoT, DgT), DgT) ds = 

c 'J c 

= \!{{P{DdT, DqT) + C{DaT, DqT, DqT), DQT)]ds. 
Substituting into (|M)) . we finally have: 

0= j {F{DqT) + {DlT,DlT)}ds, (37) 

C 

where: 

T{X):=!-R{X) + ^P{X,X)+'^C{X,X,X),x\, X eT{HTM). (38) 



Assume that J-{X) is nonnegative for any X. Then from ([37|) . it follows that 
Dg T = 0. We will show that the latter equality implies that c is harmonic. 
Indeed, taking the integral: 

= y ^ (T, DqT) ds = {{{DqT, DqT) + (T, DlT)}ds, 
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the last term in the right hand side is zero, hence {Dg^T, Dg^T) = 0. By virtue 
of the positive definiteness of g, this leads to: 

i.e., c is harmonic. We have thus proved: 

Proposition 5 //, in the Finsler space {M,F), the operator 

I 4 1 ~ 

T : T{HTM) -^ E, T{X) = ( -R{X) + -P{X, X) + -C{X, X, X), X 

takes only nonnegative values, then any closed biharmonic curve on M is a 
geodesic. 

Particular case: If g is Riemannian, then 

J^{X) = {-R{X), X) = (-7^(X, l^y)l''y, X) = -Riemg{y, n^X) ■ a, 

where Rierug is the sectional curvature of (Af , g) and a :— g{y, y)g{'K^,X, 7r*X) — 
{g{y, TT^X))^ > 0. Thus, we obtain in this case a known result, ^: any closed bi- 
harmonic curve on a Riemannian manifold with nonpositive sectional curvature, 
is harmonic. 

5.2 2-dimensional Finsler spaces 

If (Af , F) an oriented Finsler surface, the Frenet equations ([29]) reduce to: 

ei = T, Dq^ci — Kie2, -Da., 62 = — Kiei. 
and the equations of proper biharmonic curves become: 
kik']^ = 0; 



k" - kJ = (-Kii?(e2) + K?P(e2, 62) - KiKiC(e2, 62), 62) , 

that is, Ki ^ is a constant and: 

At2 = (i?(e2)-«:iP(e2,e2),e2). (39) 

Since the left hand side of the above equality is always positive, we have a 
series of immediate consequences. 

1) In Landsherg spaces, ([55)) becomes similar to the corresponding relation 
in Riemannian spaces, |12j : 

K\ = {R{e2),e2)=K{e2). 

It follows: 

Proposition 6 In 2-dimensional Landsberg spaces with nonpositive (not nec- 
essarily constant) flag curvature, any biharmonic curve is a geodesic. 

2) In locally Minkowski spaces, where R = 0, P = 0, p9l) becomes: ki = 0. 
In other words; 

Proposition 7 In 2-dimensional locally Minkowski spaces, any biharmonic 
curve is a straight line. 
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5.3 Locally Minkowski spaces 

Assume now that (M, F) is a locally Minkowski space with dim M > 3. Then, 
the biharmonic equation (22) reduces to: 

DlT + C{DaT, Dlj) = 0. (40) 

Let us choose a local chart in which gij = gij (y) and denote 

DoT=:t'S,, (41) 

the components of the tension of c in the adapted basis; with this notation, (|40p 
is read as: 



ds ds ' ^^' ds 



Lowering indices by gij and taking into account that gij.k = "^Cijk, '"' 
'i.CijkT^, this is: 



Af Dt' DTh 
ds^^'' ds ' ds 
Thus, we get the first integrals: 



ds 



= Xh, A/i = const. (42) 



So far, we have seen that two classes of locally Minkowski spaces, namely: 

- Euchdean spaces ([121); 

- 2-diniensional spaces, 

do not to admit proper biharmonic curves. Using (j42|) . we will prove that 
such a result cannot be extended to arbitrary locally Minkowskian Finsler 
spaces. 

Proposition 8 There exist locally Minkowski spaces of dimension at least three, 
which admit proper biharmonic curves. 

Proof. The space M^ equipped with the Randers-type Minkowski norm, [4], 

F{y) := y(7)^TFFTFF + hy"" (43) 

(with b G (0,1)) admits proper biharmonic curves. The detailed calculation is 
presented below, in Section 5.4.2. ■ 

5.4 Examples 

5.4.1 A 2-dimensional Randers-Numata metric 

Let us consider the unit disk M — {x eM? \ SijX^x^ < l}j equipped with the 
Randers-type metric: 

Fix, y) - \J5ijV'y^ + 5ijx\f =: a + 5^jx'v^■, (44) 

12 



and look for its proper biliarmonic curves. For this metric, we know ([1] P- 70) 
that: 

G' = ^y' (45) 

(note: our G' is as in 0, [9], i.e., half of the one in ^) and its flag curvature is: 



^i^^ify'f = 4^74- 



K - TT^iS^.yVr - —,a' (46) 



We notice that the space is projectively fiat, i.e., its geodesies are straight lines. 
The following equalities on (M, F) can be checked by direct computation: 



g,j = FF.,j + F,F, = Fa.,, + F,F, (47) 

■^ -ij — ^-ij-) -^ -ijk — ^-i'jk-) v^^j 

C^,k = \F%k = \ {F.:F.,k + F.,F.,u + F.fcF.y + FF.,,fc); (49) 

{o?)-i]k = ^ aa.ijk = -{a.ia.jk + a.ja.ik + a.ka.,j). (50) 



Along a unit speed curve c : I ^ M, s i-t- (a;*(s)) (more precisely, along its 
lift s H-S- c'{s) := (x*(s),y*(s) := i*(s)) to TM) there hold the relations: 

F[x,v) = ||r|| = l, F, = ^; F,T^^^{T,DaT)^0; (51) 

a^ 
a.,y^ = a, a.ijy^ = a.,T^ = a' + —, (52) 

where a' := da/ds = y^Sta + r^a.i and r* are as in (J4T]). The tension can be 
expressed as: 



ds +^ ^^yy - ds ^^=y - ds +^^ - ds^ f- 



In the following, we will write the equations (j39p of proper biharmonic curves 



as: 



gijT^T^ — Ki{— const. ^ 0); 

-KJ = {-R{DaT),DaT) + {P{DaT,DaT),D9T). 



(54) 



Let us examine the first equation (|54|) . Using (|47|, (fSTj) . ([531) ^-nd ([521) . this 
condition becomes: 

.,=Fa..,rr^=fa.., — — . (55) 

(y2)2 _j^ly2 (^1)2 

We have: a.n = — ^— , a. 12 = 5 — , a. 22 — — 5— and thus. 



F{yW ~ y'yY ^ id, ^ y\. 
5 ~ i'a(-— (arctan— r)) 



■ 1 1^ 
ds y^ 
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Passing to (Euclidean) polar coordinates y^ = a cos 9, y^ = a sin 61 in the plane 
(y^ ,y^), we have 9 — arctan —- + const, and the first equation (I54p can be finally 
written as: 

9 = ±Ki f{Fa)-^/^ds. (56) 

Consider now the second equation (|54p : the first term in its right hand side 
is: 

{-R(DeT),DeT) = -nl (i?(e2), 62) = -/.^F^if = 1^""- (5^) 

In order to calculate the second term {P{DqT, DqT), Dg^T) = Pj^fcrV-'r'^, 
we make use of (TlBl): 



^ijk — y ^ijk\l — V {(-^ijk,l ^ tr iL^ijk-h ^ i U^hjk ^ i jl'^ihk ^ i kl^ijhj- 



First of all, from: C tjk,i = T(-F'^;)jjfe = i^{FF,i).^jk = T(-Pai^z)jj7c, it follows, 

with the help of (gS]) and dSHl), that y^Cijk,i = 0. 

Then, taking into account that y^T\i = G'^.j, y^G\ = 2G'\ the expression 



H5|) of G" and of the fact that Cijfc is homogeneous of degree -1 in y, we obtain: 



„2 



^ijk - --^^ijk- (58) 



Substituting (H^ in (|58p . contracting with t^t^t^ and then using ([5T|). we 

h PyfcrV^r'= = 
and ([55t . we find: 



^ 2 

find: PyfcrV^r'^ = -^p.^-fcrV^r'^ = -^a.^j-^rV^r^ Finally, by 0, ([52 



(P(i.,.T,i...r),i.O.T) = M(^ + |!). 

Using the latter and (|57|) . together with F = 1 in the second equation ((54)) . this 
becomes: 

with the solution: 

a = {^ls + vf''^, 

where /i :— and i' G K are constants. Once we have calculated a, we are 

o 
able to calculate 9 (from ([55]) ). as: 



6* = ± Ki /" a-i/^ds = ±^ (^s + i^; 



3/4 , 



where 7 g R. Substituting into y^ = a cos 9, y^ = asin^ and integrating with 
respect to s, we can finally state: 
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Proposition 9 Proper biharmonic curves of{M,F) are given, in the Finslerian 
natural parametrization, by: 



x^ = J a cos 9 ds — J{fJ-s + i^)-^/^ cos(— — (^s + v) ' + 'y)ds 
x^ — J a sin 9ds — ±/(ms + i^)^/^sin(— — {fis + i/)^'* + 7)ds. 

5.4.2 A Randers-Minkowski norm on R^ 

Let us find proper biharmonic curves of the Minkowski norm: 

F{y) = Viy')^ + {y^f + {y^? + by^^:a + P (59) 

on M = R'^. That is, we will look for non-constant y = x obeying the biharmonic 
equation and F{y) ~ 1. 

Let us denote: /3 — hiy\ bi :— bS^i. The metric tensor gij is given, [4], by: 

F 

Qij = —{Sij - a..,a.j) + F.^F.j. (60) 

Also, there hold relations (|47p - ([5n)) and, along any unit speed curve c on (M, F), 
we have: 

r^ = ^, F,t' = F' = 0, a.,T' = a'; (61) 

a.,rV=F.,-rV = |, F,^ = -^, (62) 



l.„ „ ,. dr^ 

a " ' ds F 

2 



ct-ijT-' — —{SijT^ — Qfitt'), ctj— J— = C(" —; (63) 



We remark that here, nf := gijT^r^ is not necessarily constant. 
Biharmonic curves are given by P^ . i.e.: 

g^— — = Ai, Xi^ const. (64) 



Li the following, we will express these equations in terms of the Euclidean 
QCtric 6ij and of the derivatives of a. 
By dMI, dSS]), dMl), we have: 



dT^ F^ dT^ .k\ Fa'\ K? ,^^, 

For the second term in ([64]), we get, with the help of gS]), (|49]), (|6l]), (|62|) : 



C.jfe^^'^' = ^(y^- + J^a-.jfeT^'^'^) 
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and further, using (ISUl) . (|55|) : 

^--^'^ = ^P-^ + (4^ - ^)-^ - ^%-^- (66) 

Adding (|65ll and (|66|) . equations (j42]) become: 

a as 



^"'5 -1 P("')' ^"" 1 "' '^M^a -- 1 "n 


(67) 


2 "*J '2 '0 01?'' ) ^-^ '^» ' OI?"*' 

a"^ \ a^ a la It j lb 


notice some more equalities, valid along any unit speed curve: 




a.,y' = a, ^yV^ = -— = aa', %tV-' = -^ + (« ) , 


(68) 




(69) 



(where, in the third relation (|68p . we have used (|60|) ). Contracting these three 
equations (|67| in turn with y^Tr"" and (5* Afc and using (|68|) . (|69|) together with 
a.^r' = a', they take the form: 



\iy' 


= - 


nl 




XiT^ 


= K 


ik\ 




a 


dr^ 




Fa' 


ds 




a^ 



, , (F{a'f Fa" k\ k? A „, , , 

(70) 

where 7 := — tf '^*"'(Ai H ^i)Xj is ^ constant. 

2_r 

Differentiating the first equation (j70p by s, we have: A^r* = —2k\k!^, which, 
substituted into the second one yields: 

Ki = const. 



(which was found in Section 5.1. to be true along any biharmonic curve). Thus, 

we have: AiT* = =^ / 
equation (ffO)) becomes: 



dr' y^ 

we have: AiT* = =^ A^— — = 0. Further, noticing that S^^a.i — — , the third 

ds a 



1 F(aO!_Fc^ «i_^ 

a^ a2 a ^ 2a 2f'~^' ^ ' 

taking into account that F = 1 along c, this is: 

(lna)"-^ + ^+7a = 0, (72) 

which obviously admits non-constant solutions. 

Let a be a non-constant solution of (l72t . Then, substituting a = 1 — by^, 
we get a non-constant value for y^ — i^. Further, y^ — i^ and y^ = i? can be 
obtained from the relations F{y) = 1, A^y' = —k\. 
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6 Conclusions 

In the present paper, we extend the notion of biencrgy functional to curves 
on Finsler spaces and deduce the equations of biharmonic curves as the Euler- 
Lagrange equations attached to this functional; these equations turn out to 
involve, besides the flag curvature (which is the analogue of the curvature term 
in the Riemannian biharmonic equations), two specific, non-Riemannian quanti- 
ties, namely, the Cartan tensor and the Landsberg tensor of the space. Still, just 
as in Riemannian spaces, all geodesies are biharmonic curves, but the converse 
is generally not true. 

A brief analysis in the Frenet frame (which is defined similarly with the 
Riemannian case) shows that any Finslerian biharmonic curve parametrized by 
its arc length has constant geodesic curvature. 

In Section 5, we focus on the conditions under which the implication bi- 
harmonic -^ geodesic becomes true. The first case to be investigated is the 
one of smooth closed curves. In the particular case of Riemannian spaces, this 
turns out to be the known condition of nonpositive sectional curvature; still, in 
the general case, this condition is more complicated, involving the Cartan and 
Landsberg tensors. 

A special attention is paid to two-dimensional Finsler spaces (Section 5.2). 
We prove that, in dimension two, in: 1) Landsberg spaces with nonpositive flag 
curvature and 2) flat (i.e., locally Minkowski) Finsler spaces, any biharmonic 
curve is a geodesic. 

Further, in Subsection 5.3., we investigate higher dimensional locally 
Minkowski spaces. Here, we find a first integral of the biharmonic equation 
and show that, unlike Euclidean spaces, flat Finsler spaces of dimension at least 
three can admit proper unit speed biharmonic curves. 

Two examples: a 2-dimensional projectively flat space and a 3-dimensional 
locally Minkowski Finsler space, both admitting proper biharmonic curves, are 
given in Section 5.4. 
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